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Abstract 

Two components of spherically symmetric, inhomogeneous dust pen- 
etrating each other are introduced as a generalization of the well-known 
Tolman-Bondi dust solution. The field equations of this model are for- 
mulated and general properties are discussed. Special solutions with addi- 
tional symmetries — an extra Killing- or homothetic vector — and their 
matching to the corresponding Tolman-Bondi solution are investigated. 
PACS Numbers: 04.20.Jb, 04.40.Nr 

1 Introduction 

The Tolman-Bondi solution (TBS) for spherically symmetric has been 
found by Lemaitre Q 1933 and Tolman §§ 1934 as one of the first in- 
homogeneous cosmological models. Its physical properties has been widely 
investigated, e.g. Bondi 1947 ||. It has been used as a model to study the 
gravitational collapse and the formation of a black hole. In recent years 
spherically symmetric dust was studied to simulate such phenomena like 
shell-crossing or voids || . The TBS has also been taken as background 
metric for thin shells of dust to model voids and thin shells B . 

Spherically symmetric dust has been generalized in various directions. 
Mixtures of dust and radiation || or dust and a perfect fluid Q have 
been investigated. Coley and Tupper introduced a two-fluid model for 
general viscous fluids and examined various special cases. In this paper 
we consider dust shells of finite size moving with different four velocities. 
So two shells of dust can move towards each other to form a common 
region where two dust components exist simultaneously. To model voids 
spherically symmetric dust with different values for the energy densities 



has been matched by an intermediate region. Here comoving boundary 
surfaces have been considered. By introducing two dust components the 
boundaries need no longer to be comoving. Furthermore one can examine a 
central core of two dust components with a mass flux through the boundary 
that can be matched to an exterior one-component dust. 

A solution for the two-component dust (TCD) has to solve a coupled 
nonlinear system of partial differential equations (PDEs) in two indepen- 
dent variables for two functions. To find an explicit solutions further as- 
sumptions as additional symmetries are made. These special solutions still 
allow us to study different physically and geometrically interesting space 
times and lead to different topologies. 

The paper is organized as follows. In section 2 the Tolman-Bondi solu- 
tion is given. The model of two-component dust is introduced in section 3. 
The field equations are formulated in an invariant way and general prop- 
erties of the metric are discussed. Section four deals with the matching 
procedure. In the sections 5 and 6 some special solution with an extra 
Killing- or homothetic vector in the two-component dust region are ex- 
amined in detail. 



2 One component of spherically symmet- 
ric dust 

For dust, i.e. a perfect fluid without pressure, the Einstein field equations 
read 

Rab — -^Rg a b = T ab , 



with the energy momentum tensor 

T ab = pu a u b , u a u a 



1. 



0, 



(2.1) 



where p is the energy density and u a a time-like geodesic vector field. In 
a comoving spherically symmetric coordinate system the metric can be 
written in the form 



ds 2 



Y 2 (r, t)dtt 2 + R 2 (r, t)dr 2 - dt 2 



(2.2) 



where dVL 2 is the usual line element on the unit sphere. Here the solution 
of the field equations can be given explicitly, see e.g. [10]. 

For a non-comoving frame of reference, which we will use later for the 
matching procedure, the line element is best given by 



ds 2 = Y 2 dfl 2 - 



1 



Y 2 + Y a Y> 

Here Y iCL is the gradient of the sphere radius Y and 

Y = Y„u a = C U Y 



(Y a + Yu a )dx a (Y :b + Yu b )dx h - u a u b dx a dx b . 

(2.3) 



2 



the derivative of Y along u a . The free functions M and / can be expressed 
by 

M = Y(l-Y a Y' a ), (2.4) 
/ = -1 + Y 2 + Y a Y' a . (2.5) 

Their specific coordinate dependence is defined by 

C U M = M = 0, £ u f = f = Q, (2.6) 

as / and M are constant along the fluid's world lines. With the expressions 



(2.4), (2.5) and the metric ( |2.3D the field equations for one-component dust 
are fulfilled identically. 

The sign of the free function /, sometimes called the "energy function", 
leads to three different types of solutions. The hyperbolic case, / > 0, 
represents an open unbound model in contrast to the elliptic case, — 1 < 
/ < 0, which is a bound one. The parabolic case, / = 0, is marginally 
bound. The function M gives twice the mass inside a sphere with radius 
Y so it is usually called the "mass function". 



3 The model 

The energy momentum tensor for two-component dust in radial motion is 

T ab = fiv a v b + eu a u b . (3.1) 

We assume that both dust components are only coupled by their gravita- 
tional attraction and the mass density is conserved for each component, so 
they have to move on geodesies: 

v a v a = -1 v a . b v b = and u a u a = -1 u a . b u b = 0. (3.2) 

For this model one can again introduce a mass function M given by 
( |2,4|) which was first considered by Lemaitre Q and rediscovered by Misner 
for perfect fluids. M represents the active gravitational mass and two 
field equations can be reformulated with the aid of M. The rate of change 
for the mass function along each fluid line is given by 

M ta v a = eY 2 Y a (v a + u b v b u a ) (3.3) 
M >a u a = fiY 2 Y a (u a + u b v b v a ) . (3.4) 

M is only conserved for one component if the other component vanishes. 
The remaining field equations not containing the energy densities can be 
formulated in an invariant way, too. The first one requires that the tan- 
gential pressure vanishes, Geo = 0, and the second one reads 

- 2YY a , b u a v b + — u a v a = (3.5) 
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which is a consequence that u a and v a are both geodesic. 
The line element can be written as 



ds 2 = YW + 



A 2 - 1 
with 



dx a dx b (u a u b + v a v b - A(u a v b + v a u b )) , (3.6) 



u a = (0,0,m 3 , ua) v a = (0,0,1*3, Vi), 

where —A = u a v a is the scalar product of the four velocities ( \A\ > 1). 
To ensure that both four velocities are geodesic they must be gradients, 

u a = W, a v a = V a . (3.7) 

The scalars V and W are arbitrary functions of x 3 and x 4 . 

For explicit calculations two coordinate systems were quite useful. In 
the first one both velocities have only one nontrivial covariant component. 
With 

u a = (0,0,0,-1) v a = (0,0,1,0) (3.8) 
we get the line element [^] 

ds 2 = Y 2 dQ 2 + — - (dT 2 + dt 2 + 2AdTdtj . (3.9) 

In the coordinate chart ( |3.S| ) the energy momentum tensor is diagonal, 
so the field equations, G\\ = and G34 = 0, not involving the energy 
densities, are given explicitly by a a coupled nonlinear system of PDEs of 
second order: 

(A 2 + A 2 T )Y(-2A 2 - 1) + A t A T YA(A 2 + 5) 
+(A t Y t + A T Y T )A(A 2 - 1) + (Y t A T + A t Y T )(l - A 2 ) 
+ (A tt + A TT )YA(A 2 - 1) + A tT Y(l - A 4 ) 
+(Ytt + Y TT )(-(A 2 -lf) + Y tT 2A(A 2 -I) 2 = (3.10) 

A 2 {2Y t Y t + 2YY tT ) 
—A(Y 2 + YS + 2YY tt + 2YY TT + 1) + 2YY tT = 0. (3.11) 

The other non vanishing two components of the Einstein tensor determine 
the energy densities: 

G33 = H, G44 = e. (3.12) 

Another line element is of the Tolman-Bondi form, comoving to u a 
(W = -i): 

V 2 

ds 2 = Y 2 dtt 2 + v2 r ^ -dt 2 , u a = (0,0,0,-1) and v a = (0, 0, V r , V { ) 

(3.13) 

The line elements ( |3.9| ) and ( p.!3|) are related by the transformation 

T = V(r,t) t = t. (3.14) 
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In general both four velocities have got non zero shear and expansion. 
But further restrictions can be obtained by using the Raychaudhuri equa- 
tion (see [10 1 and references therein) for each velocity and considering the 
field equations. There are no solutions with positive energy densities where 
the expansion is constant along the corresponding vectors fields u a or v a 
respectively. Additionally, after a lengthy calculation one can show that 
vanishing shear of both u a and v a implies that the space time is flat 0, [D|. 

It is hopeless to find the general solution of the PDEs ( p. 10 ) and fl3.11| ). 
One way to find explicit solutions is to look for additional symmetries in 
order to reduce the number of variables. So the specific shape for the 
free functions and the radius Y of the corresponding TBS is determined 
by the matching procedure. But general properties and features of quite 
different models can be examined. Before discussing this in detail in the 
sections 5 and 6, we will consider the junction conditions between two— 
and one-component dust regions. 



4 Matching procedure 



Here we discuss the matching between solutions with the energy momen- 
tum tensors 

T ab = A* v a v b + 6 u aUb and = pu a u b 

across a time like hyper surface S. The matching procedure is based on 



the Darmois junction conditions outlined in the work of Israel [14] and 



especially for dust by Bonnor et.al. [15|. They require the continuity of 
the first and second fundamental form. E has to be tangent to one of 
the four velocities because there has to be no mass flux of the second 
dust component through the boundary surface. As a consequence of the 
junction conditions the sphere radius Y and the free functions M ( |2.4| ) 
and / Q2.5| ) of the Tolman-Bondi solution "produced" by the mass flux 
through the boundary are continuous on S. Hence M and / for one- 
component dust can be calculated completely in the two-component dust 
metric on £ because M and / are constant along the world lines of the 
corresponding dust component. Vice versa a TBS for one-component dust 
gives the initial values for the TCD on S. 

The matching of two-component dust to the vacuum is only possible 
if both velocities are parallel, i.e. in fact only for one—component dust. A 
continuous matching of the energy momentum tensor on E, that means 
one of the two energy densities vanishes, is included in the general scheme 
as a special case, too. 
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5 Solutions with an additional Killing vec- 
tor 



Special solutions with additional symmetries, i.e. Lie symmetries of the 
PDEs ( 3,1C| ) and ( 3.11| ), lead to space times which admit an additional 



Killing vector (or a homothetic vector, see section 6) with only radial-time 
components. The Killing vector £ a = (0, 0,£ 3 ,£ ) reflects a translational 
invariance in the radial-time coordinates which do not appear explicitly in 
the field equations. £ a can be space like or time like but not null because 
£ a £ a = implies e/j, < 0. Before discussing space like and time like Killing 
vectors separately in the next two subsections, we will examine common 
properties of the TCD admitting an extra Killing vector. 
In the coordinate chart ( |3,9| ) the Killing vector £ a is 

, n f 1\ , n k 2 + l-2kA . . 

r = (m,i, - s ) re. = kHA2 _ T) ■ (m) 

where the metric coefficients depend only on the reduction variable z = 
T + kt with a constant k ^ 0: 

A = A(T + kt) 

Y = Y(T + kt). (5.2) 

The field equations G\\ = ( |3.10| ) and G34 = ( p.ll| ) now read 

A 2 2k(YY" + {Y'f) 
+A(-l- (Y') 2 (l + k 2 ) -2YY"{l + k 2 )) + 2kYY" = (5.3) 
Y{A') 2 {{1 + k 2 )(-l - 2A 2 ) + kA(5 + A 2 )) 
+A'Y'{A 2 - l)(-2Jfc + A + kA) 
+YA"{{A 2 - l)(A + k 2 A) - k(A 4 - 1)) 

+Y"{(A 2 - 1) 2 (-1 - k 2 + 2kA)) = 0. (5.4) 

They have reduced to an autonomous system of ODEs where the prime 
denotes differentiation w.r.t. z. The energy densities obtained from ( |3.12| ) 
become 

u ~ E{Z) (55) 

^ " (- 1 + kA )Y 2 (A 2 -l) { j 

kE ^ (5.6) 



{-k + A)Y 2 {A 2 - 1)' 

when eliminating the second derivatives with the aid of the field equations 
(|5.3[) and ( |5.4) . Because of the common factor E(z) in the two energy 
densities ( ^.5| ) and ( |5.6| ) it is not possible that only one vanishes, i.e. the 
Tolman-Bondi solution cannot admit such a Killing vector except in the 
vacuum case. As a consequence the geometry and topology of the corre- 
sponding TBS differs from the two-dust-component. 
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When the Killing vector is orthogonal to one of the four velocities, i.e. 
k = or k = oo, one gets one-component dust which belongs to the class 



of Kantowski-Sachs because the Killing vector is space like [10|. 

In the coordinate system fl3.13| ), comoving to u a , the line element can 
be transformed in a form similar to (|5,2|): 



ds 2 = F?(z)dn 2 + ^ dr 2 -dt 2 z = r + kt (5.7) 

{kF 2 {z) + ci) 2 - 1 

with 

u a = (0,0,0,1), v a = (0,0,F 2 (z),kF 2 (z) + Cl ) and ^ = (0,0,1,-1). 

k 

(5-8) 

Because of the existence of an additional Killing vector the free func- 
tions M and / and the sphere radius Y of the corresponding solution for 
one dust component are determined completely by the Darmois junction 
conditions. So one knows the properties of the corresponding TBS by 
calculating the metric coefficients of the TCD. 

5.1 Solution with a space like Killing vector 

Solutions of two-component dust with a space like Killing vector, i.e a 
group G4 acting on a S3 (notation see fl(J), are generalizations of the 
Kantowski-Sachs class for one-component dust. Here an Euclidean "cen- 
ter" of symmetry is lacking. A radial coordinate can extend from —00 to 
00. Such models for dust and perfect fluid (where Y depends only on the 
comoving time) are also called T-models; see Ruban and references 
therein for a deeper discussion. 

Possible singularities in the energy densities might occur for Y = 0, 
which could be interpreted as a "big bang" or gravitational collapse "big 
crunch" . As an approximation of first order for a singularity at z = for 
the reduction variable one gets with the aid of the field equations 

\A\ « 1 + tX^z 2 ^, Y « cz 2/3 and e ~ — ^ fj, ' 



8ck 2 ' z 4 /3' ^ z 4/3- 

So \A\ = 1 and Y = occur simultaneously, i.e. in the neighbourhood 
of the big bang singularity both four velocities are nearly parallel. Here 
the asymptotic is the same as for one dust component with Kantowski- 
Sachs geometry. The expansion is stopped when the TCD reaches its own 
Schwarzschild sphere because the derivative Y' vanishes for Y = M when 



substituting the ansatz (5.2) in the formula (2.4) for the mass function M. 
Then the re-collapse begins and ends in the final big crunch singularity 
(see Figure 1). The T-models for dust and perfect fluid (especially the 



"uniform T-models") discussed in [16] are completely analogous; so that 



solution for itself can be interpreted as a generalization of the uniform T- 
models for two-component-dust and represents a closed universe. To get 
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a concrete solution for the TCD and the corresponding TBS the ODEs 
( |5.3| ) and (|5.4| ) have to be solved numerically. The initial values determine 
the maximal expansion, the relative velocity and the values of the energy 
densities for each dust component in a nonlinear way. 

Because of the range of Y this model can be interpreted as central core 
consisting of two-component dust related to one-component dust outside 
due to a mass flux. The mass function M is not a monotonic function (see 
Figure 2). The decreasing part indicating a mass loss for the TCD region 
is separated from the increasing part (ingoing mass flux) by the "turning 
point" Y = M. The ingoing and outgoing mass flux is not balanced 
because there is a sign change for the "energy function" / (see Figure 
2). So a part of the exterior TBS is connected to an open model while 
otherwise / lead always to a closed model (/ < 0). The properties for / 
and M can be compared to the " neck-and-two-sheets" topology with the 
same qualitative features for the TBS discussed by Hellaby [17]. 



5.2 Solution with a time like Killing vector 

For a time like Killing vector the metric is static. There exists no cor- 
responding non- vacuum solution for one-component dust. The type of 
shell-crossing singularities can occur when £ a becomes parallel to u a or v a . 
Here the corresponding energy density is diverging but the metric itself 
and the other energy density remain finite for a non-comoving coordinate 
system. Because the condition Y > M is fulfilled for time like Killing 
vectors the sphere radius Y has its maximal and minimal values when 
the shell-crossing singularities occur. To get a quantitative solution the 
ODEs Q5.3| ) and ( |5.4| ) have to be solved numerically. Numerical solutions 
show that the other energy density also reaches its extremal values at the 
shell-crossing singularities (see figure 3), e.g. 

Y\ < Y < Y 2 and m < /j, < fi 2 . 

The constants Yi, Y 2 , y,\ and pL 2 are determined by the initial values. 

Because of the finite values of Y this model can be interpreted as a thick 
shell consisting of two dust components. The shell is "produced" by two 
shells of one-component dust where the thickness of the shells is restricted 
by shell-crossing singularities. First the shells are separated then they 
crash together and penetrate through each other. So the interior of the 
"mixing zone" is static but the boundaries are varying in time. The two 
inner boundaries of the shells meet each other when the mass function M 
(figure 4) vanishes, let us say for Y = Y, because the junction conditions 
for the flat space has to be fulfilled for a moment. Before this moment the 
normal vector of the surface is orthogonal to e.g. u a then orthogonal to v a . 
Because the mass function is negative a TCD shell cannot be matched for 
Y\ < Y < Y. The shells have to re-separate because the dust components 
never become comoving. The energy function / for v a (Figure 4) changes 



8 



from positive to negative values. An ever expanding sphere is slowed down 
due to the shell crossing. That can lead to a bound state and end in an 
gravitational collapse. The other dust shell belongs to a closed model. 



6 Solutions with an additional nomoth- 
etic vector 

A homothetic vector rf represents a self-similar behaviour in the radial- 
time coordinates of the metric. Eardley jl8| examined self-similar space 
times in general. Cahill and Taub considered spherically symmetric 
self-similar space times for a perfect fluid, especially they investigated 
self-similar dust solutions. In recent years self-similar space times were 
extensively studied, for dust see e.g. [20[1 . The self-similar dust solution is 
included here as a particular solution. The main difference to the solutions 
of the previous section is that rfr] a can change sign for the same solution, 
i.e. rf can be space like or time like in different regions. 



In the coordinate system (3J}) the explicit expression for the homothetic 
vector is 

t 2 / \ 
?f = (0,0,T,t), VaV a = -^-j(l + z 2 + 2zG) (6.1) 

with the line element 

1 T 

ds 2 = t 2 F 2 (z)dQ 2 + 1 ( d T 2 + dt2 + 2G{z)dTdt) z = j- (6.2) 



The field equations (|3.1C| ) and ( 3.11 ) are thus reduced to a non-autonomous 
system of ODEs 

GG'F(z + G)(G 2 -1) 
+F{G') 2 {-\ - z 2 - bzG - 2G 2 - 2z 2 G 2 - zG 3 ) 
+G'F'{2z + G + z 2 G)(G 2 - 1) 
+G"F{(G + z 2 G)(G 2 - 1) + z(G 4 - 1)) 

+F"(G 2 -1) 2 (-1- z 2 -2zG) = (6.3) 
G + GF 2 - 2zGFF' - 2G 2 FF' 
+G{F') 2 + z 2 G{F') 2 + 2zG 2 {F') 2 
+2zFF" + 2GFF" + 2z 2 GFF" + 2zG 2 FF" = (6.4) 

which cannot be solved analytically. The energy densities can be expressed 
in the form 

■eito (6.5) 



t 2 (G 2 (z) - l)F 2 {z 



» = t 2 (G 2 (z)-l)F 2 (z) ^ (6 - 6) 
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for some functions \i\ and ei depending on the reduction variable z alone. 
The energy densities do not vanish except for the particular TBS. In the 
comoving coordinate system ( 3.1 3| ) the line element can be written in a 
form analogous to (|6.2|). 

The qualitative behaviour for the solutions of the field equations (|6.3| ) 
and (6.4) can be compared to the solutions discussed in the previous sec- 
tion. A big bang singularity with space like rf as well as a shell-crossing 
singularity with time like rf can occur both. The coordinate dependence 
of the singularities show the same asymptotic. But in contrast to the space 
like Killing vector with Kantowski-Sachs geometry a sphere can expand 
through its own Schwarzschild sphere for the closed model with big bang 
and big crunch singularity. So the mass function M is a strict monotonic 
function indicating that there always is a mass loss through the boundary 
generating the TBS outside. The phenomena with two shells of dust pene- 
trating each other and re-separating is possible, too. Here the homothetic 
vector is time like for the TCD region. 

In addition to the already known models space times admitting a ho- 
mothetic vector contain another type of solutions, presenting a compact 
ball of one-component dust inside surrounded by a shell of dust with fi- 
nite size outside. The dust shell is crashing on the dust ball and ends in 
a gravitational collapse. Here the space time is restricted by a big bang 
(big crunch) and a shell-crossing singularity where the maximal radius 
Y is reached (see Figure 5). The mass function M of the corresponding 
one-component dust is a positive monotonic function, i.e mass flux in one 
direction. The energy function / leads to an elliptic (closed) TBS w.r.t. 
to both u a or v a . 



7 Conclusions 

As a generalization of the well-known Tolman-Bondi solution for spher- 
ically symmetric dust the model of two-component dust in radial mo- 
tion was introduced. The corresponding field equations were formulated 
in an invariant way. The connection of the two-component dust to the 
Tolman-Bondi solution was discussed and the matching procedure out- 
lined. This model can describe different physical models and topologies. 
Two-component dust with a time like Killing vector can describe the cross- 
ing of thick shells of dust where the thickness of the shells is restricted by 
shell-crossing singularities. A space like Killing vector leads to a central 
core of two dust components with Kantowski-Sachs geometry and ingoing 
and outgoing mass flux. The topology of the corresponding Tolman-Bondi 
solution outside can be compared with the " neck-and-two-sheets" topol- 
ogy discussed by Hellaby. In the case with the homothetic vector the mass 
flux can go only in one direction. Additionally big bang and shell-crossing 
singularities can occur simultaneously leading to the model with a ball of 
dust surrounded by a dust shell. 
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Figure 1: Metric coefficients and energy densities for a solution of two- 
component dust with a space like Killing vector £ a = (0,0,1,-1/2) and the 
initial values Y(2) = 1, Y'(2) = -1, A(2) = 2 and A' (2) = -1. In the neighbor- 
hood of the big bang (big crunch) singularity both dust components are nearly 
comoving indicated by A = 1. 
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Figure 2: The mass function M of the exterior one-component dust is not 
monotonic because of ingoing and outgoing mass flux. The mass flux is not 
balanced because the sign change of the energy function / indicates that a part 
belongs to an open model. The energy function / belonging to the other four 
velocity is only multi plied by the factor k. 
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Figure 3: Metric coefficients and energy densities for two-component dust with 
a time like Killing vector £ a = (0, 0, 1, —2) in the coordinate chart ( |5.7| ) comoving 
to u a : with the parameter c\ = —2 and the initial values -F\(2) = 2, F[{2) = 1/5, 
F 2 (2) = —1/2 and -F 2 '(2) — 3/2. In this coordinate chart the shell-crossing 
singularity with diverging e occurs at R — where F and the other energy 
density \i has got its extremal valuesa. 
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Figure 4: Free functions for the corresponding dust component: M = belongs 
to the moment when the inner spheres of the shells meet each other. While f u , 
denoting the energy function w.r.t. u a , belongs to a closed model f v changes sign, 
i.e. an ever expanding dust shell inside is slowed down due to the shell crossing 
and can re-collapse afterwards. 
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Figure 5: Metric coefficients and energy densities with the initial values F(3) = 
1, F'(3) = 1/4, G(3) = -2 and G"(3) = -6: The scaling factor * has been set 
to unity. A big bang (big crunch) singularity occurs at Y — 0. The energy 
density e diverges because of a shell-crossing singularity when the homothetic 
vector becomes parallel to the four velocity u a . 



